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ABSTRACT 



We use the dipole expansion to provide a systematic way of including the running 
coupling into the BFKL equation. In terms of a Borel representation, we obtain 
an expression for the kernel of the BFKL equation. 
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It was first pointed out by Lipatov ]I| that the running of the coupling plays an 
important role in the BFKL equation ||. When this is taken into account the 
eigenfunctions <fii(k) of the BFKL kernel /C, given by 



d 2 k> a s K(k, k') <t>i{9) = \ <j>i(k), 



(1) 



cease to oscillate beyond a certain critical value of the transverse momentum k of 
the t-channel gluons in the "ladder" . Above this critical value the eigenfunctions 
decay exponentially and this decay can be matched to the oscillations below the 
critical value. This phase matching, together with a phase fixing for small values 
of the transverse momentum provided by the infrared dynamics of QCD, yields 
two boundary conditions which constrain the corresponding eigenvalues Aj to take 
discrete values. This results in an isolated pole for the QCD Pomeron as opposed 
to the cut obtained for fixed a s . 

The necessity of introducing the running of the coupling is at first sight surpris- 
ing since the BFKL kernel is infrared finite. This means that if the transverse 
momentum of the i-channel gluons at one end of the exchanged ladder is fixed 
by the impact factor describing its coupling to the scattering hadron, then the 
typical transverse momenta of the gluons in the subsequent sections of the ladder 
will be of the same order of magnitude. However, as explained in ref. IJ, the 
diffusion in transverse momentum along the ladder means that a wider range of 
transverse momenta contributes to the BFKL amplitude, so that the amplitude 
is indeed sensitive to the argument of the running coupling at each rung of the 
ladder. 

On the other hand, the precise mechanism used to encode the running of the 
coupling has little effect on the behaviour of the eigenfunctions at the critical 
value of transverse momentum, where the function changes from oscillation to 
exponential decay. In other words, it makes little difference whether the argument 
of the coupling at a particular rung is set equal to the transverse momentum above 
or below that rung. In terms of the integral equation this means that one has the 
freedom to set the argument of the coupling equal to the external (unintegrated) 
value, k, rather than the integrated one k', and write eq.(|l|) as 



The reason for this is that although far enough along the ladder we expect to have 
to consider transverse momenta whose magnitude differs substantially from the 
value set by the impact factors, the transverse momenta of two adjacent sections 
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of the ladder are indeed of the same order. As far as (renormalisation group 
improved) perturbation theory is concerned, therefore, a discussion about which 
transverse momentum should be used in the running of the coupling will only 
affect the solutions to the BFKL equation at the subleading logarithm level. 

However, if one wishes to consider in more detail the infrared contributions to 
BFKL and run the coupling down to the region of the Landau pole, where small 
changes in the argument lead to substantial changes in the value of the coupling, 
one no longer has the above freedom and a study of the precise prescription for 
running the coupling becomes pertinent. Such a precise prescription is necessary 
for an analysis of the renormalon structure of the BFKL amplitude aimed at 
identifying non-perturbative power corrections to the amplitude which are likely 
to play a crucial role in the reconciliation of the perturbative QCD Pomeron with 
the phenomenological "soft" Pomeron. 

In a recent publication || Levin has conjectured that the correct prescription for 
the coupling is 

a s (k' 2 )a a ((k-k') 2 ) 

His analysis, however, was based on the idea of considering the running of the cou- 
pling at one rung of the ladder only, whilst keeping the coupling fixed throughout 
the rest of the ladder. A systematic treatment requires a prescription which can 
be extended to the entire ladder. In this letter we propose just such a systematic 
treatment in order to arrive at a well-defined algorithm for the introduction of 
the running of the coupling. 

Our starting point is not the gluon ladder approach of BFKL but rather the 
dipole expansion of Mueller and Nikolaev and Zakharov ||. The difference 
between the two approaches lies in the description of the same high energy (or 
small x) physical process. In the BFKL approach the process occurs via the ex- 
change of a ladder of reggeized gluons in the t-channel. The leading logarithmic 
corrections in the energy are generated by the increasing number of rungs along 
the gluon ladder. The exchanged gluons, though, cannot be interpreted in a prob- 
abilistic fashion as emanating from a specific hadron involved in the process. The 
gluon ladder does not have a straightforward parton model interpretation with 
an identifiable hard scattering section. In the dipole approach the process occurs 
via the sequential emission of colour dipoles in the s-channel by the participating 
hadrons. The leading logarithmic corrections in the energy are generated radia- 
tively by the increasing number of emitted dipoles. This dipole cascade can be 
interpreted in the appropriate gauge (light cone gauge) as a component of the 
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Figure 1: Graphs for the emission of a soft gluon with transverse momentum k\ 
(impact parameter y[) from a dipole with momentum p, in which the quark has 
transverse momentum q (impact parameter x). 

hadronic wavefunction. The hard subprocess is encoded in the dipole cross section 
that describes the scattering of two dipoles for onium-onium elastic scattering, 
or the absorption by the nucleon target of a dipole emitted by the virtual photon 
for small x D.I.S. In impact parameter space, where transverse separations x and 
x' replace the transverse momenta k and k', the factorization of the dipole cross 
section from the squared wave function for dipole emission takes a particularly 
simple form. 

The equivalence of the two approaches for fixed a s and for inclusive observables 
such as onium-onium total cross section and F2 in small x D.I.S. has been shown 
in refs. [f| and @ respectively. For a systematic inclusion of the running coupling 
we find the dipole approach more suitable exactly because of its radiative nature. 
Our method is the following. The emission of a daughter dipole by a parent dipole 
in momentum space occurs via the radiation of a soft gluon with transverse mo- 
mentum k\ as shown in fig.l. We take the argument of the running coupling 

^2 

in this case to be equal to k\ , in accordance with the treatment of soft gluon 
radiation in the case of timelike cascades ||. We note that the relevant kine- 
matic region for the momentum of the soft gluon is the Glauber region where its 

-> 2 

virtuality is determined by its transverse momentum, i.e. k\ « —k\ . To trans- 
late this running into impact parameter space representation we require that the 
Borel transformed probability for emission of a soft gluon with impact parameter 
y[ from the parent dipole with size x should be equal to the expression obtained 
by first Borel and Fourier transforming the two graphs of fig. 1 and then squaring 
their sum. Since this procedure correctly reproduces the singularity structure in 
the Borel plane, it provides a prescription for the coupling that allows us to run 
down towards the Landau pole. It turns out that the argument of the running of 
the coupling cannot be expressed simply as a function of x and y[ only, but the 
expression we obtain nevertheless has a convenient integral representation. 

To this end we consider the emission of a soft gluon with momentum k%, colour a, 
polarization e and longitudinal momentum fraction z±, from a fast moving dipole 
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with total momentum in which the quark and antiquark have momenta q^ 
and (p — k\ — qY and m q is the quark mass, 



q^ = {pz 



m q + q 2 

pz 



Q). 



Leading logarithmic contributions in the energy are obtained in the kinematic 
region of strongly ordered rapidities z\ <C z. The amplitude for such a process 
with running coupling as described above is given by M 



^ a (q,h,z, Zl ) = -2g(h)T a \^°\h,z) - ^((h + q),z) 



ki ■ e 



>2 i 



(3) 



where ip^ is the wavefunction for the dipole without gluon emission and T a is 
a colour generator in the adjoint representation. We now introduce the Borel 
transformed amplitude ip^ a defined by 



^ a (q,k 1 ,z,z 1 )= / d6^ (1)a (g,fci,z,^i,6)e- fe/as(M2) . 
J o 



(4) 



Note that the Borel dual variables are b and a s (M 2 ), where M 2 is an overall 
external momentum scale to be determined once the dipoles are considered in the 
context a specific process. For D.I.S. at small x, for instance, M 2 is the spacelike 
virtuality Q 2 of the photon 7*, which fluctuates into the dipole cascade before 
the interaction with the nucleon target. The Borel transform of the wavefunction 
in eq.(P|) is 



^ a {q,k 1 ,z,z 1 ,b) = -2g(h,b)T a \^°\h,z) - ^\{k x + $,z 



ki ■ e 



ki 



• 2 ' 



(5) 



where 



—* z 

g{h ,b) 



Vb \M 2 



2\ -b/3o 



(6) 



and (3q = ( l/47r) ( (11/3) N c — (2/3) Nf). The Borel transformed g is determined by 
the requirement that the "square" of g(k\ , b) in the sense of Borel convolution 
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—* z 

with itself results in the Borel transform of 47ra s (ki ), i.e. 



^a s (h ,6) = ( 



^2 ^2 

,6) ® <pi ,6) 



" db'gih^b'Mh 2 , (6-6')). 



(7) 



Taking the Fourier transform of eq.(f|) we obtain the wavefunction as a function 
of impact parameters x, yi conjugate to q, k\ respectively 



^ a (x,y-[,z, Zl ,b) 



1 2i r°° 
-T a ^°\x,z) / drJ 

V6 7T JO 



1 T 



X 



yi \ e • vi , / M - yi ) \ 



+ 



e • (x - l/i, 
(x - yi) 2 



where Ji is the first order Bessel function. The Borel transform, $(x, z, 6), of the 
total probability for the emission of a gluon from a dipole separated by impact 
parameter x in which one of the fermions carries a fraction of longitudinal mo- 
mentum z, is obtained by taking the "square" of the amplitude given by eq.(H). 
By "square" we mean convolution of ■j/n 1 )" in 6 with itself (as defined in eq.(|7])) 
and integration over the phase space of the outgoing gluon. This gives 



4> (1) (x>,6) 



C F r z dzi 



dVi / drdr , J 1 (r)J 1 (rO 



x 



x 



TC J Z\ J JO 

]_ \ub/3 f \ \ (l-w) 6 A) (M 2 ) &/3 ° 



dcu 



1/1 \ I 

Co> 2 (l — Cl>J 2 



7" 



/2 



yl 2 (f - y!) 2 



;(yi 2 ) 6ft (f - yl) 2 + ((x - yl) 2 ) 6/ Vi 2 

+2yi ■ (x - y[)(yx 2 r bPo ((S - $(»)(£, z ) 



(9) 



where $ < -°- ) is the probability for the dipole without gluon emission. The integra- 
tion over the variables r, r' can be readily performed. However if we choose not 
to do this then each term in the integrand can be identified as the Borel transform 
of a running coupling with a suitable argument. Thus we define 

ri 2 H = (j& 2 ) w ((if-j&) a )< 1 - w > 



X 2 (t,t'u) = (r 2 ) £J (r ,2 )( 1 -^, 
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and take the inverse Borel transform of eq.(|9|) to obtain the following expression 
for the probability for gluon emission with running coupling 



& x \x,z) = -4ln(0)/dVi/ dTdr'h{r)h{r') 

TT J JO Jo UJ2[\ 

[a s (\ 2 /fi 2 )x 2 + (a s (A 2 M 2 ) - a s (A 2 /ri 2 )) (x - y[) 2 
+ [a s {\ 2 / {x - y[) 2 ) - a s (X 2 /fl 2 )^j yl 2 ] &°\x,z). (10) 



y x 2 (x -yi) 2 



In the limit where the coupling is fixed the integrations over r, r' , u) give a factor 
of 7i and the first term on the R.H.S. of eq.([3I]) coincides with the expression 
obtained in the usual treatment of dipole radiation (see, for example eq.(10) 
of ref. 0|) whereas the other two terms vanish. For running couplings we are 
unambiguously led to the prescription of eq . (|10|) , namely that the kernel has 
three terms each of which has a different running coupling and for each of these 
terms the argument of the running coupling is not expressed simply in terms of 
the impact parameters x and y, but also in terms of three further parameters, 
r, t', u, the emission probability being a weighted integral over these parameters. 
It is worth noting that the weight function peaks when A(r, r',u) is equal to 
unity, so that although one has to integrate over all possible arguments for the 
running coupling, the integral is highly peaked in the region where the running 
depends only on y 2 or (x — y[) 2 or some mean of the two. We note that had we 
taken the Fourier transform of the wavefunction ip^ a in eq.^) and squared it to 
construct as in eq. (|T0"D , we would have ended up with a product of the form 
g(r/y 2 )g(T' / (x — y) 2 )- It is the use of the Borel representation that enables us 
to express this product in terms of running a s of a single argument. This then is 
a consistent and systematic way to introduce the running of the coupling, valid 
to all orders in perturbation theory. 

So far we have dealt with the wavefunction squared $ without considering its 
coupling to some hard scattering subprocess. Such a coupling was considered 
in ref. || where, through the use of time ordered perturbation theory, it was 
shown that final state interactions cancel and are therefore not included in the 
dipole kernel. Consequently, final state interactions will not affect the method 
we suggest for the running of the coupling. 

Returning to the Borel transformed $W, in eq.(|9|), the integrations over y{, r, r', 
to may be performed to yield 

^n {XiZib) = ._ Hz) __^_j (ii) 
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which has the usual infrared renormalon poles for positive integer values of bfto 
[§]. The poles at b = and for negative integer values of bflo are ultraviolet 
renormalon poles and are eliminated if a short distance cutoff is introduced for 
the integration over y[. Unitarity tells us that the one gluon exchange virtual 
contribution is —^^(x,z,b). 

If we square the amplitude for single gluon emission calculated in momentum 
space, take the Borel transform of the running coupling and then take a Fourier 
transform in the transverse momentum q of the fermion only we get upon inte- 
grating over all phase-space for the outgoing gluon 



¥ 1] (x,z,b) = %ln(z) 



7T- 



d 2 h 
h 2 




2 — g* fc l'^ _ Q~ tKi 



fcl -X 



x,z) 



(12) 



which also yields eq.([[l]) when the integral over k\ is performed. This is not 
surprising but it serves as a check and confirms that one is at liberty to commute 
Fourier and Borel transforms. 

In exact analogy with the treatment of ref.[[| the process can be iterated to 
account for the emission of any number of soft gluons. To this end we define 
the generating functional Z (x, 0, z,j(y' , z'),bj to be the quantity whose n th 

functional derivative with respect to the source, j(y',z'), generates the Borel 
transform of the probability for the emission of n gluons. Due to the expo- 
nentiation of soft gluon radiation, virtual corrections can be accounted for by 
the factor exp ( — z, ||. Then the generating functional Z obeys the 
Bethe-Salpeter equation 



Z (x, } z,j,b 

c F r*dz' 



2C F . r (-&#,) m*x 

exp | — — m(z 



2^2 \ b /V 



71 



+ - 



7T" 



i , 

— exp | m(z/z 



T(l + b(3 )\ 4 

T(l + bf3 ){ 4 ) 

1 du / 1 \ wb ^° ( 1 \ ( 1_aJ ) b '' 3 o 



f d?y' f°° drdr'hirW) f , , 

J -JO JO UJ2 (1 — uj)l \T 

M 2b/3 



T 



12 



y (x-y y 



(yy^(x-y') 2 + ((x-y') 2 )^y' 2 

+2y' ■ (x - y' )(y' 2 r b ?°(x- y' )*)W* 

g> Z (y", x, z'j, bj ® Z (y\ 0, z',j, b) j(y', z'), 



(13) 
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where the exponentials of functions of b are to be understood in the sense of their 
Taylor expansion with products replaced by convolutions, i.e. 



exp (/(&)) = 5(b) + /(&) + -/(&) <g> /(&) + 



The n <?i functional derivative of eq. flHf ) with respect to the source j is equivalent 
to the n th iteration of the integral equation 



o ~ o f°° ~ ~ 

z—$(x, z, b) = z— $ (0) (x, z)S(b) + / dx'Kix, x' , b) <g> z, b), (14) 
oz oz Jo 

where JC is the Borel transformed dipole kernel and a change of variables has been 
affected from the two dimensional vector y' to its modulus y' and the modulus, 
x', of x' = (x — y'). The Jacobian for this change of variables can be written as 



d 2 y' = 2n y 1 x'dy'dx' I ndnJofax) Jo(kx') Jo(m/)> (15) 
J o 

where Jo are Bessel functions of zeroth order. Extracting /C from eq.(|i~3|) and 
integrating over y' and k, we find 

C A 1 (M 2 x' 2 \ bPo /"i du T(l - (1 - w)&/3 ) r(l - ubf3 Q ) 



n 2 x'\ 4 J Jo r(l + (l-w)6A))r(l+a;6A)) 



2 \ ^fti^ 1 / 2 ' 

x < |% F 1-%,1-%;1;^ 



/ r 2 _ „/2\ / 2 \ <"6A> / „2 ~ 

(2 \ w&ft) / 2 1 

F^ ,-^ ;l;^]^, (16) 

where x< = min(a;, x'), x > = max(x, x') and F is the Gauss hypergeometric 
function. We have replaced Cp by \C^ which is the colour factor occurring in 
the BFKL equation and agrees with Cp in the leading 1/N C approximation. The 



8 



first term of eq.(|l^) describes the virtual corrections and it is proportional to 
8(x — x') due to the completeness of the Bessel functions. 

The eigenfunctions, 4>(x, b) of this kernel are also functions of b so that the action 
of the kernel is also to be understood in the sense of a convolution, i.e. the 
eigenvalue equation is 

poo pb „ _ 

/ dx' / db'JC(x,x', b')<j>i{x', (&-&')) = Aif&Wz, b). (17) 

JO JO 

The eigenvalues and eigenfunctions of this kernel must be found in order to inves- 
tigate the analyticity structure of the Borel transform of the BFKL amplitude. 
Clearly eq.flT7D does not lend itself to a simple analytic solution. The eigenfunc- 
tions of the kernel will have to be examined by numerical methods in order to 
locate their singularities. The location of these singularities will give us informa- 
tion about the non-perturbative contributions to the BFKL amplitude. We shall 
report on the results of such an analysis in a forthcoming publication. 
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